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Abstract We present a relatively simple time domain method for determining the 
bandpass response of a system by injecting a nanosecond pulse and capturing the 
system voltage output. A pulse of sub-nanosecond duration contains all frequency 
components with nearly constant amplitude up to 1 GHz. Hence, this method can 
accurately determine the system bandpass response to a broadband signal. In a novel 
variation on this impulse response method, a train of pulses is coherently accumulated 
providing precision calibration with a simple system. The basic concept is demon¬ 
strated using a pulse generator-accumulator setup realised in a Bedlam board which 
is a high speed digital signal processing unit. The same system was used at the Parkes 
radio telescope between 2-13 October 2013 and we demonstrate its powerful diag¬ 
nostic capability. We also present some initial test data from this experiment. 

Keywords Astronomical instrumentation, methods and techniques — Methods: 
bandpass calibration, pulse calibration, nanosecond pulse generator 


1 Introduction 

Precise calibration of the bandpass of a radio telescope is critical for a number of as¬ 
tronomical observations, including searches for the spectral signature of the cosmo¬ 
logical evolution of neutral hydrogen EH and for baryon acoustic oscillations ffl. 
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The bandpass of an instrument is most commonly calibrated in the frequency do¬ 
main, by measuring the spectrum of an injected calibration signal, but it can also be 
calibrated in the time domain, by injecting a time-varying signal, such as a pulse, and 
recording the pulse profile after it has passed through the instrument. The recorded 
pulse profile provides an indication of the impulse response of the instrument, and its 
Fourier transform — divided by that of the original pulse — represents the complex 
spectral response. 

In this work, we extend this impulse response technique by generating a regu¬ 
lar series of pulses with timing locked to our sampling rate, allowing us to coherently 
sum the profiles of the individual pulses. Compared to using individual pulses or cali¬ 
bration with an incoherent broadband noise source, this allows us to rapidly extract an 
accurate measure of the impulse response of the instrument, while keeping the ampli¬ 
tude of individual pulses sufficiently low to avoid non-linear behaviour. The required 
processing is minimal compared to frequency-domain methods utilising spectrome¬ 
ters, which must calculate the spectrum in real time. Implementation is also simpler 
than for a similar technique used in telecommunications in which a known pseudo¬ 
random sequence is recorded instead. The generated pulses are well-measured and 
substantially shorter than the inverse bandwidth of our system, constituting a well- 
understood flat-spectrum calibration signal. The total power in individual pulses is 
negligible relative to the system noise, so that the pulse injection system could be run 
continuously without degrading the system performance. Our approach constitutes a 
simple technique for accurately calibrating the complex spectral response of a radio 
telescope. 

In section [2] we describe the hardware used in the implementation of our tech¬ 
nique, consisting of a fast pulse generator and a Bedlam board 0 used for digital 
signal-processing. In section [3] we describe a benchtop implementation of our tech¬ 
nique in a system similar to the SARAS spectrometer (6), in which the signal propa¬ 
gation paths and other non-idealities are well understood, to demonstrate its function¬ 
ality. In section]?] we show some results from the implementation of our technique on 
an active astronomical instrument, the Parkes radio telescope. Finally, we conclude 
briefly in section [5] 


2 System description 

The basic experimental setup for system bandpass calibration is shown in Fig. |T] We 
inject a nanosecond pulse to a system being calibrated. System responses to several 
pulses are coherently added in an accumulator. The system components are described 
below. 


2.1 Nanosecond pulse generator 

The pulse generator we used for this work and the output pulse it produces are shown 
in Fig. [2] This pulse is injected into the system being measured and the system output 
is sampled. 
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Fig. 1 Schematic of the pulse accumulator setup. 



Fig. 2 Nanosecond pulse generator (Left). A typical output pulse voltage seen in an oscilloscope (Right). 
The horizontal ticks are 0.1 ns. Sampling rate of the oscilloscope is 20GS/sec and bandwidth is 6 GHz. 


Since radio telescope receivers have high gain and sensitivity only a low voltage 
pulse is required, as it is desirable to keep the input voltage excursion of the pulse 
below the input noise power level, to ensure the system retains a linear response to 
the pulse input. The pulse is generated with a custom-built circuit based on a high 
speed positive emitter coupled logic (PECL) D-type flip-flop with reset. The pulse 
generator board works in two modes: Asynchronous triggering from an on board 
time base or externally triggered from a supplied event. The generated pulse has the 
following properties: 

- Output pulse width : 350 ps (50% width). 

- Output rise/fall time : 240 ps. 

- Pulse amplitude : 400 mV unipolar pulse into 50 Q. 
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The low intrinsic jitter of 0.2 ps for the PECL D-type flip-flop allows synchronous 
integration of the received pulses to be performed to high precision as this is neg¬ 
ligible compared to the sample period of the following Analog to Digital Converter 
(ADC), which is approximately 1 ns. The narrow pulse width and rise/fall time en¬ 
sure the pulse contains spectral content up to a few GHz, which was sufficient for our 
purposes. The power requirement for the pulse board is 3.6 V at 100 mA. 

As our injected pulse is not a perfect delta function, having non-zero width and 
some asymmetry, strictly it is necessary to deconvolve the injected pulse from the 
received pulse profile in order to achieve absolute phase and amplitude calibration. 
However, as the pulse width is substantially narrower than the inverse of the 300- 
350 MHz system bandwidth in sections [3] and [4] it is close to a perfect delta function 
at these timescales, and so we omit this step. 


2.2 Bedlam Board 

In our experiments, the trigger signal is sourced from a custom general purpose high 
speed digital signal processing board, known as Bedlam. The Bedlam board consists 
of 8 RF input channels feeding two 8-bit quad-channel ADCs (EV8AQ160). The 
digital sample data from each of two channels of the ADCs are directed to one of four 
DSP oriented FPGAs (XC5VSX95T). A final logic oriented FPGA (XC5VLX30T) 
provides dual fast Ethernet and general purpose I/O interface from the board. 

The pulse generator output is fed into the system being calibrated and the system 
output is captured by the Bedlam board (Fig. [TJ. In order to improve the signal-to- 
noise ratio using the coherent summation of pulses in the pulse train, the sample 
voltages are accumulated inside the Bedlam FPGAs. The accumulator buffer length 
is programmable with a buffer length of 8192 samples used in the experiments de¬ 
scribed here. A sampling rate of 1.024 GHz is used such that each buffer is 8 ^sec 
long. The pulse generator is then triggered at a rate of 125 kHz so that one pulse 
is produced every 8 /I sec. The trigger is synchronised with the pulse generator so 
the pulse voltages will be added coherently. The Bedlam board is programmed to 
dump the accumulated buffer using an Ethernet interface to a general purpose com¬ 
puter, when the required accumulation count is reached. At the pulse repetition rate of 
125 kHz, and duration of 1 second, 1.25 x 10 5 pulses are accumulated resulting in a 
350:1 improvement in the signal-to-noise ratio of the captured pulse. Eonger integra¬ 
tion periods can be used if desired. Initial tests showed some repetitive weak additive 
signals caused by sampling clock leakage into the ADC. Since these are clearly syn¬ 
chronous with the sampling clock they were removed by differencing the two halves 
of the output buffer with the pulse present in only one half, but the clock contamina¬ 
tion equal in both halves. The Bedlam board functional schematic and accumulator 
block diagram are shown in Figs. [3] and [T] respectively. 

Fig. [4] shows a typical measured pulse, accumulated over 1 second, after differ¬ 
encing the two halves of the accumulation buffer. Only the first hundred samples are 
shown to show the output pulse profile. This is the pulse, fed directly to the accu¬ 
mulator, without the system being measured inserted in between the pulse generator 
and the accumulator. The Bedlam board has a voltage inversion at the input due to a 
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balun that converts the single-ended to differential, so the pulse amplitude is negative 
at the output of the accumulator. The polarity of the pulse, however, does not have a 
significant effect on the bandpass calibration. 



Fig. 3 Bedlam board schematic. 
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Fig. 4 Accumulator output after 1-second integration of the 1 MHz pulse train. 
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3 Experimental setup for bandpass measurements 

To validate our experimental setup for bandpass calibration, we performed a benchtop 
test of a system similar to the SARAS spectrometer 0. This system contained a 
low-pass filter, and compared signal paths with and without a T connector with an 
open load, to replicate the bandpass ripple observed in SARAS measurements. The 
block diagram of the system under test is shown in Fig. [5] with the pulse generator 
connected at its input and the accumulator connected to its output. Two channels of 
the Bedlam board are used with identical accumulator setup and common triggers. 
After injection, the short pulse is band limited by a lowpass filter with a cut off of 
350 MHz. Filter output is split into two halves by a power splitter and one is fed 
directly to an accumulator via channel B of the Bedlam board output which is used 
as a reference. The other half is fed to channel A via a T connector to mimic internal 
reflections that are typically present in receivers, and which must be calibrated. The 
third input of the T is connected to a 10 dB attenuator via a 5 m long cable and the 
output of the attenuator is connected to an open load. A part of the incident wave 
at the input of the T reached channel A of the Bedlam board directly whereas the 
other part is reflected from the attenuator output and appears as a delayed pulse at the 
accumulator input in channel A. 

Fig. [6] shows the direct and the delayed pulses at the output of the accumulator 
in channel A after 1 second integration. Band limiting the signal to 350 MHz results 
in ringing of the pulse in the time domain. Fig. [6] shows the reference pulse in chan¬ 
nel B and the pulse received in channel A. Since the pulse power is divided between 
the direct and the reflected pulse, the direct pulse amplitude is smaller in channel A 
compared to the same in channel B. 



Fig. 5 Experimental setup for pulse calibration. The power combiner, low-pass filter (LPF), splitter and 
T connector together constitute the system under test. 
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Fig. 6 Comparison of the channel A and channel B output after 1 second integration of the 1 MHz pulse 
train. Both channels contain the primary pulse which are identical and coincide in time. The pulse in 
channel A, however, goes through the T junction. A part of the pulse energy is lost due to the loss in the 
T junction and also a part of it is transmitted through the open end of the T which returns to channel A 
after reflection as a delayed pulse at a delay of 175 ns. Hence, the primary pulse amplitude in channel A is 
smaller compared to that in channel B. 


The true spectrum of the band limited pulse along channel B with the spectrum of 
the pulse along channel A are shown in Fig. [ 7 ] Introduction of the T and long cable in 
channel A have two obvious effects. In the time domain the delayed attenuated pulse 
is clearly seen with a delay of 120 ns and this caused the fast (16 MHz) ripple in the 
channel A voltage amplitude spectrum. The voltage amplitude spectrum also shows 
a 320 MHz ripple which correspond to a 6 ns delay and is caused by reflection of the 
direct pulse between the Bedlam board input and the T output. The 6 ns delay is less 
conspicuous in the time domain as it is confused with the pulse oscillation caused by 
the sharp cut-off in the low pass filter. 

As our results are consistent with the properties of the system under test, we 
consider this to be a successful benchtop validation of our bandpass calibration setup. 
Assuming the input pulse to be sufficiently narrow (see section ing , Fig.[7]then shows 
the amplitude of the system bandpass. 


4 Pulse calibration tests with the Parkes radio telescope 

An opportunity arose to test the system described in section [2] when three Bedlam 
boards were used at the Parkes radio telescope between 2-13 October (2013), in an 
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Fig. 7 Band limited output spectra. Y-axis is the Fourier transform of the voltage time samples. 


experiment intended to detect pulsed radio emission from atmospheric particle cas¬ 
cades initiated by cosmic rays or gamma rays. The Bedlam boards were originally 
designed for the LUNASKA experiment (4), in which an alternate mode was used 
to detect and store individual pulses coming from space. This experiment required 
high-precision timing calibration between the beams of the Parkes 21 cm multibeam 
receiver Q in order to discriminate between radio-frequency interference, which typ¬ 
ically appears in multiple beams simultaneously, and atmospheric particle cascades, 
which are expected to move from beam to beam with a timescale of nanoseconds. 
This calibration was achieved with the same pulse generator system described in 
section [2] so we were able to carry out coherent pulse accumulation with the same 
hardware. 

The gross physical experimental setup is shown in Fig. [8] The signal from the 
pulse generator was transmitted from the apex of the parabolic dish, received at the 
focus, and ultimately returned to the Bedlam board, similar to Fig.[l] The system un¬ 
der test in this case comprises the large-scale structure of the telescope, the multibeam 
feed, and the remainder of the receiver signal path, including low-noise amplifiers, 
downconversion and attenuation. Transmitting the pulse from the apex of the dish 
was convenient and flexible for this experiment, but may differ in some respects from 
observations of astronomical sources, which are in the far field of the telescope. 

The pulse generator was set to produce pulses at a cadence of 125 kHz. Fig. [9] 
shows the accumulated pulse profile after one second of integration. The direct pulse 
is most prominent, but there are eight additional pulses also visible, with successively 
smaller amplitudes, with a spacing of ~ 180 ns. These result from the partial reflec¬ 
tion of the transmitted pulse, first from the base of the focus cabin, then from the 
vertex of the dish. These points are separated by 26 m, with a round-trip light travel 
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Fig. 8 Experimental setup for pulse calibration tests with the Parkes radio telescope. The pulse is trans¬ 
mitted in the band 0.7-4 GHz from the (red) log-periodic dipole antenna (LPDA) at the vertex of the dish 
surface, and received by the 21 cm multibeam receiver at the focus with a narrower bandwidth (BW) of 
300 MHz centred on 1.38 GHz. Some fraction of the pulse is reflected (yellow) from the focus cabin, and 
is responsible for the reflected pulses described in the text. The received signal is transmitted by cable 
(dashed) from the focus to a Bedlam board in the control tower. (Schematic by J.M. Sarkissian.) 


time matching the observed pulse spacing. In the frequency domain, these reflections 
create a ripple in the bandpass of the system, which is a well-known problem seen in 
most radio telescopes; for Parkes, the ripple has a period of 5.7 MHz IZI71 . match¬ 
ing the inverse of the pulse spacing. The higher-order reflections (i.e. the second and 
subsequent reflected pulses) will result in harmonics which introduce sharper ripples 
in the system bandpass. A more detailed analysis of the reflected pulses reveals a 
complicated pattern of reflected polarized radiation which results from the structures 
around the telescope focus cabin. 


5 Conclusion 

We have explored a time domain technique for calibration of wideband spectrome¬ 
ters, which is particularly well suited for applications in radio astronomy. This tech¬ 
nique involves injecting a series of pulses into a system input, and coherently accumu¬ 
lating the system output, with no further need for complicated real-time processing: 
the amplitude and phase of the system bandpass response can be determined from 
the Fourier transform of this accumulated output. The key advantages of this tech¬ 
nique are that the total power contained in the individual pulses is negligible relative 
to the system noise, allowing it to run continuously without degrading system per¬ 
formance; and that it allows both the amplitude and phase of the bandpass response 
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Fig. 9 Pulse profile recorded at Parkes using the pulse generator-accumulator system. The spacing between 
the successively weaker pulses indicates that they result from repeated reflection of the pulse between the 
focus cabin and the vertex of the dish (see Fig. [8}. 


to be measured with a system implemented on simple hardware. We have tested an 
implementation of this technique for bandpass amplitude calibration in a benchtop 
experiment, and with the 21 cm multibeam receiver system on the Parkes radio tele¬ 
scope. 

The science goals of various proposed radio-astronomy experiments will require 
precise calibration of the system bandpass, which further development of our tech¬ 
nique can in principle achieve. For example, detection of the global signature of the 
epoch of reionization in the background radio spectrum requires a calibration accu¬ 
racy of 1:10 5 . Detection of baryon acoustic oscillations, spectral features imprinted 
on the redshifted 21 cm background at z ~ 1, requires a calibration accuracy of 1:10 6 . 
Detection of the recombination of primordial plasma to form neutral hydrogen, which 
causes a broad spectral feature in the spectrum of the cosmic microwave background 
around 1-20 GHz, requires a calibration accuracy better than 1:10 9 . While this last 
case is beyond the time resolution possible with the current implementation of our 
technique, developments in fast sampling and signal processing may enable this in 
future experiments. 

Our technique also has a potential application in the calibration of ultra-wideband 
receiver systems, which are currently being planned and constructed for several tele¬ 
scopes (e.g. 0). These receivers use a single feed with a bandwidth of, typically, 
several GHz, which is then split into several overlapping sub-bands before being sam¬ 
pled and digitised, to improve robustness against strong radio-frequency interference 
and to reduce the sampling-rate requirements. After digitisation, these sub-bands are 
combined in order to restore the complete, continguous frequency band, but this re- 
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quires calibration of the phase and amplitude differences between these sub-bands. 
The coherency of our pulse calibration system allows this to be done, as the phase and 
amplitude of the pulse components must be continuous across the sub-band bound¬ 
aries. 
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